THE WEAK AND STRONG CLOSURES 
OF SOBOLEV HOMEOMORPHISMS ARE THE SAME 



TADEUSZ IWANIEC AND JANI ONNINEN 

Abstract. Let X and ¥ be bounded multiply connected Lipschitz domains 
in R 2 . We consider the class J^(X, Y) of homeomorphisms h: X °" to > Y in 
the Sobolev space W 1,v {%., R 2 ) . We prove that the weak and strong closures 
of Jtf p (X, Y) , 2 ^ p < oo , arc equal. The importance of this result to the 
existence theory in the calculus of variations and anticipated applications to 
nonlinear elasticity are captured by Theorem 1.5. 



1. Introduction 

The primary theme of this paper is about homeomorphisms h : X "-^^ Y between 
planar domains in the Sobolev space # /1 ' p (X, Y) . Throughout we use the notation 
=^,(X, Y) for the class of such mappings and reserve the symbols Jf p (X,Y) and 
Jf p (X, Y) for the closures of Jf p (K, Y) in strong and weak topology of W 1,P (X, Y) , 
respectively. 

Theorem 1.1. LetX andY be l-connected Lipschitz domains. If £ > 2 and p ^ 2, 
then 

(1.1) JT P (X,Y) = j%(X,Y). 

Equality (1.1) also holds for I — 1 if p > 2 , but fails when p = 2 . 

For simply connected domains the equality J^2(X,Y) = ^2(X,Y) also holds 
under suitable normalizations of the mappings h : X Y . This case is fully 
resolved in §4.3. 

It should be noted that the weak closure Jf^,(X, Y) is none other than the set of 
all weak limits of homeomorphisms h : X Y in the Sobolev space !# /1 ' p (X, Y) . 
This fact is not obvious but follows from Theorem 1.4, see Remark 2.7. 

The interest in Sobolev homeomorphisms comes naturally from Geometric Func- 
tion Theory (GFT), Calculus of Variations (CV), and from Nonlinear Elasticity 
(NE) where prospective applications to elastic plates and thin films gain additional 
motivation. In GFT we seek, as a generalization of the celebrated Riemann Confor- 
mal Mapping Problem, a homcomorphism h G ^^,(X, Y) with smallest p -harmonic 
type energy integrals. It is certainly unrealistic to expect that the infimum energy 
within the class ^(XY) will always be attained; injectivity is often lost when 
passing to the weak limit of the minimizing sequence. The best example is the 
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collapsing phenomenon in the minimization of the Dirichlct integral for mappings 
between circular annuli [4] . It is new, even for the basic case of the Dirichlet energy, 
that 

(1.2) inf ■ [ [ \Dh(x)\ 2 dx = inf ff \Dh(x)\ 2 dx . 

At this point let us stress that the Dirichlet energy integrals are quintessential in 
the study of harmonic mappings [20, 23, 29, 37, 38, 39]. In GFT all known examples 
of mappings with smallest energy appear to be strong limits of homcomorphisms. 
This observation is now a fact confirmed by Theorem 1.1. 

In the proof of Theorem 1.1 we first consider mappings h which are weak limits 
of homeomorphisms h j : X Y in the Sobolev space # /1 ' p (X,Y). They need 
not converge strongly. Nevertheless, for the proof of Theorem 1.1, it seems natural 
to try suitable corrections in the sequence hj to gain strong convergence. This 
did not work. The strong approximation will be achieved by modifying the limit 
map h instead. The first observation in our approach is that the weak limit of 
homeomorphisms hj extends continuously to a monotone map h : X Y between 
the closures. Monotonicity, the concept of C.B. Morrey [47], simply means that the 
prcimage hr 1 (C) of a continuum in C € Y is a continuum in X. The main ingredient 
in the proof of Theorem 1.1 is a series of step by step replacements of h by piecewise 
p-harmonic diffeomorphisms. In making such replacements we heavily rely on the 
weakly converging sequence of homeomorphisms hj : X '^^>- Y . WE believe that 
once the monotonicity of h is established or assumed the sequence {hj}JL 1 should 
play no role in the construction of strong approximation of h . This rises even more 
general question. 

Question 1.2. Let X and Y be planar domains. If ft,: X Y is a monotone 
Sobolev mapping in #' 1 ' P (X, Y), can h be approximated strongly in # /1 ' p (X, Y) by 
homeomorphisms hj : X Y? 

Uniform approximation of monotone mappings with homeomorphisms is of great 
interest in topology. We refer to Rado [52] and Youngs [59, 60] for the earliest con- 
tributions and to [45, 53, 57] for further information. In the Sobolev setting the 
approximation problem is at the very heart of GFT, and is important in mathe- 
matical models of elastic deformations. Indeed, Question 1.2 is closely related to 
the Ball-Evans question: 

Question 1.3. [10, 12] If h e #' 1 ^ P (X, Y) is invcrtible, can h be approximated in 
# /1 ' p (X, Y) by piecewise affine invertible mappings? 

This question is of great interest in the study of neohookean energy functionals [8, 
14, 18, 24, 54]. The Ball-Evans problem for planar bi-Sobolcv mappings which 
are smooth except for a finite number of points has been resolved in [46]. For 
analogous questions in the planar Holder continuous setting see [15], and for bi- 
Lipschitz homcomorphisms see [21]. The difficulties in saving injectivity in the 
process of piece- wise linear (equivalcntly smooth if n = 2,3) approximation are 
already recognized in [10, 51]. In [30] we showed that every homeomorphism h: X — > 
Y between planar open sets that belongs to the Sobolev space #' 1,P (X,Y), 1 < 
p < oo, can be approximated in the norm topology of this space by ^"^-smooth 
diffeomorphisms, See also [31] for the approximation up to the boundary. 
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In this paper we actually prove stronger statement than Theorem 1.1. Namely, 
we obtain an approximation by diffeomorphisms which have the same boundary 
values as h . Precisely, we have 

Theorem 1.4. Let X and Y be (-connected Lipschitz domains, 2 I < oo, and let 
hj : X Y be homeomorphisms converging weakly in W 1,P (X, Y), 2 ^ p < oo, to 
a mapping h e W 1,P (X,Y). Then there exists a sequence of ( ta" x -diffeomorphisms 

h* :X°-^Y, h*eh+ W^ P (X, Y) 

converging strongly in W 1,P (X,Y) to h. This also holds for simply connected Lips- 
chitz domains if p > 2. 

Theorem 1.4 has applications to a study of traction free problems in NE [6, 7]. 
Suppose we are given a well defined energy functional 

(1.3) £ x [h] = J J E(x, h, Dh) dx , for Sobolev mappings h € W 1 ' P {X, Y) . 

Here the stored energy function E:XxYxR 2x2 — »■ R + satisfies the usual 
Caratheodory regularity conditions and coercivity E(x,h,£) >p |£| p , p > 2. How- 
ever, the key hypothesis on E is the Morrey's quasiconvexity [48], which just 
amounts to the lower semicontinuity: 

(1.4) £ x [h] < liminf£x[/ij], whenever hj h weakly in W 1 ' P (X,Y) . 

We shall actually restrict (1.3) to the subclass Jf p {X,Y) C #' 1 ^(X, Y) of weak 
limits of homeomorphisms. Recall our standing assumption (Lipschitz regular- 
ity) on the domains X and Y which is sufficient to have a continuous extension 
h : X Y of h e Jf? p (X,Y) together with a uniform bound of the modulus 
of continuity by means of the energy. This is immediate from Lemma 2.9 and 
Lemma 2.10. Now it makes sense to speak of the partial boundary condition for 
h € J#p(X,Y) 

(1.5) h(x) = h Q {x) , for x e X C dX 

where h , referred to as a boundary data, is a Sobolev homeomorphism in Jf p (X, Y) 
and X is a closed subset of dX. Condition (1.5) is void when the set X C dX is 
empty This latter case, being energy-minimal analogue of the Riemann Mapping 
Problem, became of great interest in GFT. Another important borderline case, 
which might be of interest in NE, occurs when X = dX . There is, however, a 
marked difference between our setting and the classical boundary value problems 
in NE. Here we find the minimizers in the class J^ p (X,Y) , which is in the closet 
proximity to homeomorphisms. The Lagrange-Euler equations are no longer avail- 
able [9, 11, 55]; they have to be replaced by the inner-variational equations; also 
known as energy-momentum or equilibrium equations, etc [19, 50, 56]. Various 
results in this direction were obtained in [13, 20, 32, 46, 58]. 

In the general case of (1.5) we say that h is traction free on dX\X , see [6, 7]. We 
consider the family si — ^e(/i ) C Jf? p (X,Y) of homeomorphisms which satisfy 

(1.5) to inquire into its infimum energy. 

(1.6) £(*0 =mf{£ x [h] ; hesJ x (h )} 

Let si and si denote the closures of si in strong and weak topology of the Sobolev 
space #' 1 ' P (X,Y), respectively. The essence of Theorem 1.4 is now captured by 
the following result 
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Theorem 1.5. Let X and Y be Lipschitz domains. Under the assumption on the 
energy integral, specified above for equation (1.3), there exists f £ srf C Jf p (K,Y) 
such that 

(1.7) // E(x,f,Df) dx =E{sf) 



Here one particular case, when X and Y are simply connected and the exponent 
p = 2 , requires that the set X C <9X contains at least three distinct points, see 
Section §^.5. 

In particular, letting X — , we obtain 

(1.8) £x[f] = inf (( E(x,h,Dh) dx = inf [ [ E(x,h,Dh) dx 

Some comments about the proof of Theorem 1.4 merit mentioning here. First 
we conveniently reduce the problem to the case of domains bounded by circles 
(Schottky domains). Surprisingly this takes some efforts to justify such a reduction; 
namely, continuity (in the Sobolev norm) of the composition operator with a bi- 
Lipschitz transformation has to be secured. The main idea of the proof, like in [30] , 
is to partition X into small cells U C X and replace h in U by a p -harmonic 
diffeomorphism. However, the novelty lies in the replacements over the boundary 
cells, see §3.3. We explore topological properties of monotone mappings to see the 
geometry of boundary cells; they are certainly simply connected domains. Then, 
with the assistance of the weakly converging sequence {hj} , we construct a p - 
harmonic diffeomorphism having the same values on the boundary of a cell as h . 
Here, for p > 2 , a p -harmonic variant of the classical Hurwitz theorem might be 
of independent interest. 

At this point it may be worth mentioning our definition of a p-harmonic mapping. 
Let h = u + iv be a complex- valued function in the Sobolev space #^ 1,P (X, C) , 1 < 
p < oo . We consider the gradient map Vh = (Vu, Vv) : X '-^^ R 2 x R 2 and its 
«S? p -energy. 

(1.9) S x [h]= J J \Vh(x)\ p dx = J J (\Vu(x)\ p + \Vv(x)\ p ) dx 

The motivation for choosing exactly this form of the ^f p -energy functional comes 
from the Lagrange- Euler system which in this case consists of two independent 
scalar p-harmonic equations 

(1.10) div|Vw| p - 2 Vu = and div| Vv\ p - 2 Vv = 

Definition 1.6. A complex-valued function h = u+iv of Sobolev class W lo ' c p (n, C) , 
1 < p < oo , is called a p -harmonic map if it satisfies the uncoupled system of 
equations (1.10). 

This notion is different from what can be found in the literature. We take 
advantage of Definition 1.6 by not actually appealing to the coupled p -harmonic 
system div \Dh\ p - 2 Dh = , see §3.2. 

2. Notation and Prerequisites 

2.1. Domains. Throughout this text X and Y will be, unless otherwise stated, 
planar Lipschitz domains of finite connectivity 1 ^ I < oo . Thus each boundary 
<9X and dY consists of I disjoint closed Jordan curves which are locally graphs of 
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Lipschitz functions. We reserve the following notation of the ordered I -tuples of 
boundary components, 



We shall often distinguish the outer boundaries of X and Y. These are the 
boundaries of the unbounded components of C \ X and C \ Y , respectively. An 
^-connected domain X is said to be Lipschitz regular if each boundary component 
Xi,X2,...,Xf is a closed Lipschitz curve, locally (upon rotation) the graph of a 
Lipschitz function. Let us point out that any ^-connected Lipschitz domain can 
be transformed via a bi-Lipschitz mapping : C C onto a domain whose 
boundaries are circles (Schottky's domain). However, it is not so apparent whether 
there is such ^ whose composition with Sobolev functions represents a continuous 
(nonlinear) operator. Section 2.5 is devoted to constructing this type of bi-Lipschitz 
transformations. 

2.2. The class Jff(K, Y) of homeomorphisms. We shall consider orientation 
preserving homeomorphisms h: X Y . In general, a homeomorphism h: X 
Y may not extend continuously to the closures. Nevertheless, each boundary com- 
ponent of Y is a limit set under h of one and only one boundary component of 



2.2.1. Boundary Correspondence. The boundary components of X and Y can be 
so numbered that the limit set of 3L V is T v . Let us assume this numbering and 
record it as 

(2.2) h: X hence for the inverse map h : T v ~-> 3£„ , v = 1,2, ...,£ 

This notation precisely means that for every e > there exists S > such that 
dist[/i(x), T„] < e, whenever dist[x, £„] < S. 

Definition 2.1. Given two ^-connected domains X and Y , each with the specific 
order of the boundary components as in (2.1). We denote by Jf? (X, Y) the class 
of orientation preserving homeomorphisms h: X Y satisfying (2.2). 

Remark 2.2. Without losing any argument in this paper, one could have restricted 
our definition of Jf(X, Y) to homeomorphisms in a designated homotopy class. 

2.3. Monotone mappings. The notion of monotone mappings has originated in 
the work of C.B. Morrey [47]. 

Definition 2.3 (monotonicity) . A continuous mapping h: X Y between topo- 
logical spaces is said to be monotone if for each y e Y the set h^ 1 (y) C X is a 
continuum; that is, compact and connected. 

Since we will be dealing with continuous monotone mappings between compact 
subsets of the extended complex plane C = C U {0} , it is appropriate to recall a 
theorem of Whyburn [45] . 

Theorem 2.4. A continuous mapping /: X Y between compact metric spaces 
is monotone if and only i// _1 (C) C X is connected for each connected set C C Y. 



(2.1) 




for the components of <9X 
for the components of dY 



X 



We also recall a theorem of Kuratowski and Lacher [41, 42]. 
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Theorem 2.5. Let X and Y be compact Hausdorff spaces, Y being locally con- 
nected. Suppose we are given a sequence of monotone mappings fa : X Y 
(such are homeomorphisms) converging uniformly to a mapping f: X — > Y, then 
f: X Y is monotone. 

For further reading about monotone mappings we refer to [45, 53, 57]. 

Proposition 2.6. Let f: C C, / : C ' 2£i;i > C 7 &e a continuous monotone map. If 
IcC is a continuum that disconnects C into two components then so is 

Proof. One needs only observe that there are two types of points in C\/ _1 (K): the 
ones that can be path-connected with oo (they make the unbounded component) 
and the ones that cannot be path-connected with oo. These latter points are 
mapped by / onto the bounded component of C \ K, thus form the remaining 
(bounded) component of C \ / _1 (K). □ 

2.4. The classes Jf p (X,Y) , i^(X,Y) and ifp(X,Y). The following classes of 
mappings of finite energy will be used throughout. 

• Jf p (X,Y) = if(X,Y)n:T 1 ' p (X,Y). 

• JZp(X,Y) is the closure of Jf p (X,Y) in strong topology of W 1>P (X,C)- 

• Jf p (X,Y) is the closure of Jf p (K,Y) in weak topology of /T 1 ' P (X,C). 

• ^ P (X) is the Royden p- algebra of uniformly continuous functions on X 
having finite p-harmonic energy. The norm in ^ P (X) is given by 

IHU(x) = ™PlM*)l + ( JJ\VHx)\ p tey 

• M P (K) is the completion of ^^(X) in the above norm. 

Remark 2.7. Some reflections concerning the terms weak and strong closures are 
in order. If if is a subset of a Banach space W , then its weak closure if is the 
smallest subset of W that contains if and is closed under weak convergence. The 
strong closure if , on the other hand, consists exactly of strong limits of sequences 
in if . Thus if C if . However, in general, when if is unbounded then the set 
of all weak limits of sequences in if need not be weakly closed. Consequently, if 
is larger than the set of all weak limits. However, by virtue of Theorem 1.4, it is not 
difficult to see that ifp(X, Y) consists exactly of all weak limits of homeomorphisms 
h : X Y in the Sobolev space #' 1 ' P (X, Y) . It is for this reason that we do not 
introduce a separate notation for the class of weak limits of homeomorphisms. The 
interested reader may consult our subsequent paper [36] for further applications. 

Remark 2.8. Functions in Royden p-algebra extend continuously up to the bound- 
ary, because of uniform continuity. Thus we identify the Royden p -algebra as 

(2.3) ^ P (X) = Sf (X) n W^{X) 

This is a commutative Banach algebra with multiplication and addition defined 
pointwise. Note too, that if X is sufficiently regular (such are the Lipschitz do- 
mains) then / e ^ P (X) belongs to ^(X) if and only if f(x) = on <9X . It 
is also clear that if a sequence {fj}jLi of functions fj e ,^ P (X) converges to / 
uniformly on X and weakly in ^^(X.C) then / e M P (X) . 
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2.4.1. The class Jf p (K, Y),p > 2. Since X and Y are Lipschitz domains, their 
boundaries are neighborhood retracts so we infer from [35] that 

Lemma 2.9. Every h € Jf p (K,Y) extends as a continuous monotone map h: X°-^ 
Y. Moreover, the boundary map h: dX <9Y is also monotone. 

Lemma 2.10. We have the uniform bounds of the modulus of continuity. For 
fteJf 2 (X,¥) and 1^2, 

(2.4) \h( Xl ) - h(x 2 )\ 2 < g 2 (X,Y) _ |vfe|2 ^ a;ijX2G - 

log 



|xi-x 2 | y 

for any bounded Lipschitz domain X C C and /i G ' P (X, C), p > 2, we have 

(2.5) \h{ Xl ) - h{x 2 )\ p ^ C p {X)\ Xl - x 2 \ p ~ 2 J J \Vh\ p , x u x 2 eX. 

In view of (2.4) and (2.5) a sequence {hj} C Jf p (X,Y) that is converging 
weakly to h actually converges uniformly on X. Furthermore, such homeomor- 
phisms hj : X Y are monotone so the boundary mappings hj : <9X dY are 
monotone as well. Therefore, the weak limit h: X Y is also monotone and 
h : <9X dY is monotone. It is worth noting that the inequalities (2.5) and (2.4) 
remain valid for h £ J0p(K,Y). In particular, we have a continuous imbedding. 

(2.6) i^(X,Y) C^ P (X) 

2.5. Transition to circular domains. Let Y be a bounded domain in R 2 and 
f:Y°^Y'a bi-Lipschitz deformation so its inverse F -1 : Y' °-^> Y is Lipschitz 
as well. The induced composition map 

(2.7) F t : W lj, (X,Y) -> W^ P {X,Y'), F t (u) = F o u 
and its inverse 

(2.8) Fj" 1 : ^•"(X.Y') -)■ ^ llP (X,Y), F," 1 ^) = F" 1 o w 

are well defined nonlinear bounded operators, where X can be any domain in R 2 . 
It will be advantageous to transform Y into a domain whose boundary components 
are circles, so-called circular or Schottky domain. We will require the induced 
composition operators Fj and F^ 1 to be not only bounded but also continuous. 
Unfortunately, this is not always the case , see [26]. Our goal is to construct, a 
special bi-Lipschitz transformation of Y into a circular domain. 

Proposition 2.11. Given any bounded ^connected Lipschitz domain Y C R 2 , 
there exists a bi-Lipschitz map F : R 2 R 2 which takes Y onto a circular domain. 
The induced composition operators Fj : #^ 1,P (X, R 2 ) — > W 1,P (H, R 2 ) and the inverse 
Ff 1 : #^P(X,R 2 ) -> #' 1 'P(X,M 2 ) are continuous for all 1 < p < oo. 

Let us comment on some related results. First consider the Sobolev space 
#' 1:P (X,M) of real-valued functions. If /: R — > R is Lipschitz the composition 
f o u with u e /^ /1 ' P (X,M) represents a Sobolev function whose gradient can be 
defined by the rule 

^. ( f'(u(x)) Vu if / is differentiable at u(x) 

I otherwise. 
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The point is that / is diffcrcntiable everywhere except for a set £cKof zero linear 
measure. The preimage u~ 1 (E) C X may have a positive Lebesgue measure. But 
in this case Vw vanishes on this set. We refer the reader to a paper by Marcus 
and Mizel [44] in which they show that in fact the induced operator 

(2.9) /,: #' 1 'P(X,M)^^ 1 ^(X,M), /,(tt) = /o« 

is continuous. Now suppose F: Y Y' is a Lipschitz map that is ^-smooth 
except for a set E C Y of isolated points. Then the induced composition map 

(2.10) F t : W hp (X,Y) -> #' 1 ' P (X,Y'), F i (h)=Foh 
with h € W 1 ' P (X, Y) is still well defined. The differential is given by 



D(Foh) 



DF(h{x)) Dh(x) , whenever h(x) & E 
, otherwise. 



The preimage h 1 (E) may have positive measure but Dh(x) vanishes on this set 
anyway, ft is not difficult to see that the operator (2.10) is continuous. 

Proof of Proposition 2.11. We first approximate Y by polygonal domains. Let T v 
be one of the boundaries of Y. Consider a partition of T„ into closed Jordan 
arcs A\A 2 , A 2 A%, . . . , A n -iA n , A n A n+ \, A n+ \ = A\ , defined by a sequence of 
consecutive points A\, A2, . . . , A n G Y„ . Associated with such a partition is a 
polygonal chain with vertices at Ai, . . . , A n ; that is, a piecewise linear curve P v 
that consists of the line segments A\A 2l A 2 A 3 , . . . , A n _iA n , A n A\ connecting 

n — 

1 he endpoints of the arcs. Although the union T„ = (J AiAi + i is a closed Jordan 

i=l 

curve the polygonal chain may have points of self-intersection. However, since T„ 
is Lipschitz regular there is e > small enough so that if 

(2.11) dist(^i, A 2 ) < e, . . . , dist(A„, A{) ^ e 



then the only consecutive segments A\A 2l ^2^3, A n _\A n , A n A\ intersect, 
of course at their common endpoint. Under the condition (2.11) we obtain the 
boundary of a simply connected domain, called polygon. With e > sufficiently 
small we may, and do, ensure that each arc AiA 2 , . . . , A n Ai , becomes upon a 
rotation a graph of a Lipschitz function; points slightly above the graph lie in Y 
and points slightly below the graph lie in M 2 \ Y . Choose and fix one of the arcs 
, say the graph of a Lipschitz function y = f(x) , a ^ x ^ b . Here, for 
some small positive 8 , we have 

Q+ = {(x, y):a^x^b, f{x) < y sC /(*) + 5 } C Y 

SI- = {(x, y) : a < x < 6, f(x) -S^y^ f(x) }cR 2 \Y. 

With such 8 fixed we may, if needed, further partition the arc AiA i+ i into a finite 
number of consecutive subarcs 

A~A~l 1 = A~B[ U BJh U . . . U B^V^i 



so that each straight line segment AjBi, P>iP> 2 , B^Aj + i lies strictly below 
the graph of the function y = f(x)+8 and strictly above the graph of the function 
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y = f(x) — S. However, to simplify the writing we assume without loss of generality 
that the arc AjAj + i already enjoys this property. Thus we have the following region 

ft = 0+ U fT = {(x, y): a^x^b, f(x) - 5 sC y sC f(x) + 6} 

and two cross cuts with endpoints Aj = (a, /(a)) and A/+i = (&,/(&)). One cross 
cut is the graph of y = f(x), a < x ^ b, and the other cross cut is a straight line 
segment AjAj + i 7 which we shall view as graph over the interval [a, b] of a linear 
function <p = <p(x), 

(x-a)f(b) + (b-x)f(a) 



V = 



for all x e 



b — a 

Note that f(x) — S < ip(x) < f(x) + S for a < x < b. It will be convenient to extend 
/ as a Lipschitz function in the entire real line R, simply by setting f(x) = ip(x) 
outside the interval [a, b}. We are now in a position to define a bi-Lipschitz mapping 
$: M 2 -» R 2 (associated with the arc AjAj + i C dY) as follows, 

<&(x, y) — (x, y') for — oo < x < oo, — oo < y < oo 



where 



y = y + 



ipjx) - f{x) 
25 



(\y - f(x) -S\- 2\y - f(x)\ + \y- f(x) + S\) 



y, 

f(x)+S-ip(x) , 

s y 

<p(x)+6-f(x) , 

s y ~t~ 



f(x)+S 
S 

/(*)- 



y > fix) + s 

[tp(x) - f{x)] , fix) sc y < f{x) + 6 
[fix) - <pix)] , fix) -S^y^ fix) 
y < fix) - s 



We see from this later formula that if x is fixed the function y' = y\x, y) is piece- 
wise linear and strictly increasing in y. Thus $ is a bi-Lipschitz map. Moreover, 
$ix,y) — ix,y) for all ix,y) e R 2 \ ft, which is immediate when y ^ fix) + S 
and y ^ fix) — S . On the other hand if x (a, 6) then <p(a;) = /(x), so y' = y 
from the first formula. Also note that <& takes the arc AjAj +1 into the line segment 
connecting Aj and A,+i; indeed, for y = fix) we have <E>(x, fix)) = ix,ipix)). 
An important feature of this particular map $ : ] 
continuous operator $j = ^ 1,P (X,R 2 ) ^ ^^(X, 
hkiz) = Uk + Ufc(z) converging strongly in #' 1 ' P (X,! 
have 



: 2 -> R 2 is that it induces a 
I 2 ). Indeed, given a sequence 
x 2 ) to hix) = u(x) + ivix) , we 



Mfc, «fe + 

— ► u, w + 



2(5 

ipju) - /(m) 
25 



(K - fiu k ) -6\- 2\v k - f{u k ) \ + \v k - fiu k ) + 5\) 
i\v - fiu) - (5| - 2\v - fiu)\ + \v- fin) + 5\) 



The final step in the construction of the map F : R 2 — > R 2 consists of dividing the 
entire boundary dY into a finite number of sufficiently small arcs Ti, T 2 , . . . , C 
dY. To each arc I\ there corresponds a bi-Lipschitz map $ l : R 2 ]R 2 , i — 
1, . . . , N , which takes 1^ into a line segment connecting the endpoints of Tj . Note 
that $ J is the identity map on all remaining arcs. Therefore, the composition 

F := $ 1 o • • • o $ w : R 2 ^> R 2 
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is a bi-Lipschitz map which takes Y onto a polygonal domain and the induced 
operator 

F B = $Jo---o*^: W 1 - P {X,R 2 ) ^ W 1 ' V {X,R 2 ) 

is continuous. The same applies to the inverse map -F -1 ; it induces a continuous 
operator F^ 1 as well. Lastly, we compose F with a bi-Lipschitz map G : R 2 — > R 2 
which is a -diffeomorphism outside the corners of Y' . This results in a map 
G o F: R 2 — >• R 2 that takes Y into a circular domain. 

Remark 2.12. We may choose the circular domain such that the center of the outer 
boundary lies outside the domain. This of course requires that I > 2 . The centers 
of the remaining boundary circles certainly lie outside the domain. This additional 
requirement will later help us to perform the reflections about the boundary circles. 
The case 1=1 will be treated differently in §4. 

Remark 2.13. An analogous transition from the domain X into a circular domain 
X' , via a bi-Lipschitz map T : C "-^^ C , T : X X' , presents no difficulty. 
The induced composition operator 

T": ^^(X'.Y) >T 1 ' P (X,Y) T s u = uoT 
and its inverse are always continuous, 1 < p < oo. 

□ 

From now on, unless otherwise stated, both X and Y are 
circular domains containing no centers of the boundary circles. 

2.6. Extension to C. For convenience and also for easy references, we shall extend 
h: X Y to a monotone mapping of C = C U {oo} onto itself. Here, with the 
usual convention concerning the point oo, the one point compactification C = C U 
{oo} will be identified (topologically) with S 2 . Recall that we are working under the 
assumption that the boundary components X\,. . . ,X( C <9X and Ti, . . . , Tg C <9Y 
are circles. To each circle Xi = {z: \z — Zi\ — ri} there corresponds the reflection 
map 

<Pi-. C->C , tpi(z) = Z Z% r 2 + Zi , i = l,...,£. 

\Z Zi\ 

Denote by Xj = (pi(X) the reflected domains. These domains together with X 
furnish a circular domain of connectivity 1(1 — 1) 

X+ = X U Xi U • • • U Xg l) X . 
The boundary components of X + are denoted by 

X] = tpi(Xj), i^j. 
Similarly to X, we consider the boundary circles Tj of Y and define the corre- 
sponding reflections ipi : C — >• C. We obtain mutually disjoint domains Y, = ipi(Y), 
i = 1, . . . ,£, which together with Y furnish a circular domain 

Y+ = Y U Yi U • • • U Y e 1) Y . 

The boundary of Y + consists of £(£ — 1) circles T\ = ^(Tj), i ^ j. Now, let 
h : X Y be any continuous mapping that takes each Xi C dX into the corre- 
sponding boundary component Ti C <9Y . We also assume that in this correspon- 
dence h takes the outer boundary of X into the outer boundary of Y. 
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2.6.1. The extension h : X + Y+. There is a natural way to extend h contin- 
uously to a map of X + onto Y + . The notation h : X + Y + will be reserved for 
such an extension of h : X — > Y . By definition, the map h : X + Y + takes the 
reflection of X into the corresponding reflection of Y 

h = ipi o h o ipi : Xj — > Y, for i = 1, 2, . . . , I . 

It should be noted that if h : X — > Y is monotone and continuous, then so is 
the mapping h : X + Y + . Since in our applications the boundary mappings 
h: Xi °-^>- Tj arc monotone it follows that h: X\ , i 7^ j, arc monotone 

as well. Furthermore, whenever h: X Y is a homeomorphism, the extended 
mappings h : Xi Yj, i = 1, . . . ,£ , are homeomorphisms. 

We will explore the Royden p-algebra. Note that if h e ^(X, Y) n #' 1 ' P (X, Y) = 
^p(X,Y) and h: X t =*2> T;,fori = 1, then ft e ^(X+, Y+)n# /1 ^(X+, Y+) = 

^" P (X + , Y + ). We have a uniform bound of the p-harmonic energy 

// |Vfc| p <C p (X,Y) // IV/iP*. 

2.6.2. Weak continuity of the Jacobian determinant. If a sequence of mappings 
fk € #^ 1,P (X,C), p > 2, converges to / weakly in W 1,P (X) then the Jacobian 
determinants J(x,fk) converge to J(x,f) weakly in Jz? x (X) . We write it as 

(2.12) [[ <p(x)J(x,f k )dx — > ff ip(x)J(x,f)dx, for every (p e JSf°°(X) 
J Jx J Jx 

However, if p = 2 , this property is invalid for several reasons. The best example to 

illustrate is a sequence of Mobius self-homeomorphisms of the unit disk converging 

to a constant map, see §4 formula (4.1). Note that in this example all the Jacobians 

arc nonnegative, yet the integrals (2.12) fail to converge. The situation is quite 

different if the mappings in question admit W 1,2 - extension beyond X with a 

nonnegative Jacobians. Precisely, we have 

Lemma 2.14. Suppose X is compactly contained in a domain X + ; such is our 
Schottky domain X and its extension X + by reflections. Let a sequence of map- 
pings fk ■ X + '-^^ C with nonnegative Jacobians be bounded in W • (X + ,C) and 
converge to f weakly in W 1 ' 2 (X,C) . Then (2.12) holds. 

As a particular case, if t/ C X is a measurable set we obtain, upon setting 
<P = Xu ! tnc following variant of (2.12) 

(2.13) JJ^J(xJ k )dx — ► JJ^J(xJ)dx 

In other words the averages of J(x, fk) over any measurable set U converge to the 
average of J(x, f) . The proof is not straightforward, but relies on the well known 
££ logJz? - integrability of nonnegative Jacobians [49]. The interested reader is 
referred to [16, 34] for further reading about Jacobians. 

2.6.3. An auxiliary extension h : C C . Proceeding in this direction, we shall 
further extend h : X + Y + to the entire Riemann sphere C , again denoted by 
h : C C . Let us point out that we will need this extension only to provide easy 
references to the theory of monotone mappings on § 2 ; for example, Theorems 2.4 
and 2.5. Let X^ and Y? , i ^ j , denote the disks enclosed by the circles X\ and 
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Tj , respectively. There is one exception to this notation; namely, the case when 
the circles %[ C <9X + and Tj C <9Y + are the outer boundaries. In this case X^ 
and stand for the complements of the disks. We then extend each boundary 
map h : X\ Tj , i ^ j , in a radial fashion into X^ and Y^ , and continue to 
write h: Yj, i ^ j . As an illustration, the radial extension into the disks 

X? = {re 10 : sC 6 < 2tt} and t{ = {Re l9 : < 6 < 2tt} takes the form 

h: X{^Y{ , h(pe 10 ) = P -h(re ie ) 

for ^ p ^ r, and for r < p ^ oo in case of the complements of the disks. 

There will be no need for Sobolev regularity of h outside X+. Of course, con- 
tinuity of h : C — >■ C is understood with respect to the choardal metric induced by 
the stereographic projection of C onto § 2 . 

Observe that a uniform convergence of continuous mappings hj : X Y, 
hj : Xi — > Tj, i — !,...,£, j = 1,2,... to a mapping h: X Y yields uni- 
form convergence of hj : C 2 ^ 2 > C to h: C C. If, in addition, all the mappings 
hj are monotone then so is the limit mapping h . 

2.7. The p -harmonic Dirichlet problem. There are two common settings of 
the Dirichlet problem. The classical one, with a continuous boundary data, com- 
bines Perron method and Winer's criterion of regular points on the boundary of a 
domain. In the variational approach, on the other hand, one seeks to minimize the 
energy integral over the class of functions in u + Wo' p {^l) , where u G W 1 ' p (Vt) 
is viewed as the boundary data. Even when these two different settings are well 
defined the question whether they lead to the same solution involves a delicate 
analysis of the boundary of the domain. Strangely, in the widely spread theory 
of the Dirichlet integral, explicit statements concerning simply connected domains 
appear to be rare in the literature. The equivalence of these two settings is vital in 
our approach. 

Theorem 2.15. Let U C R 2 be a bounded simply connected domain and u e 
^(9)f]W 1 ' p {Vl), 1 < p < oo. Then there exists unique function u e ^(Ti)nW h P(n) 
that is p-harmonic in and equals u on dfl . Furthermore, u € u + Wo' p (Q) and 

£ p [u] < £ p [u] 

equality occurs if and only if u = u. 

Proof. We shall not give all details for the proof. Nevertheless, it is worth remarking 
(because it is not obvious) that the variational solution u extends continuously up 
to the boundary. This is because each boundary point of a planar simply connected 
domain is a regular point for the p-Laplace operator A p [25, p. 418]. See [28, 6.16] 
for the full discussion of boundary regularity and relevant notion of capacities. In 
particular, we refer the reader to [27, Lemma 5.3], see also [40, Lemma 4.1] and [43, 
Lemma 2] for a capacity estimate that applies to simply connected domains. □ 

2.8. Univalent p-harmonic extension. The celebrated Rado-Kneser-Choquet 
theorem asserts that a harmonic function h : Cl — > C in a Jordan domain, which 
extends continuously as a homeomorphism of dQ onto a closed convex curve L C 
C , is a -diffeomorphism of fi onto the bounded component of C \ T . In fact 
our proof of Theorem 1.4 in case p = 2 relics heavily on the theory of harmonic 
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mappings. An excellent reference is [23] . Similar arguments will apply to the general 
case of p > 2 . However, new ingredients (the p -harmonic variant of Hurwitz 
theorem) will be needed. Meanwhile, let us call upon the following variant of the 
Rado-Kneser-Choquet theorem from the work of Alcssandrini and Sigalotti [2], see 
Theorem 5.1 therein. 

Theorem 2.16. Suppose that C C is a simply connected Jordan domain, G a 
bounded convex domain and h € ^(fi) H a p -harmonic map that takes 

dQ, homeomorphically onto dG . Then h : G is a < ^°° -diffeomorphism. 



3. Proof of Theorem 1.4 for multiply connected domains 

Some additional prerequisites, more specific to our proof, are in order. We 
present them in the following subsections. 

3.1. Squares and cells. Let for a moment h: C C be a general continuous 
monotone map, h: C — > C, h(oo) = oo. We shall work with an open square QcC 
and its preimage = h~ 1 (Q). The set f2 is called a cell, which is surely a domain. 
Its complement C \ 0, being equal to ft. -1 (C \ Q), is also connected. In view of 
unicoherence of C, the boundary <9f2 = fin(C\fl) is connected. Thus f2 is a simply 
connected domain. As for the preimage h~ 1 (dQ) a caution is required. Although 
it is straightforward that dh _1 (Q) C /t -1 (dQ) the set C := h _1 (0Q) can be signifi- 
cantly larger than dh~ 1 (Q). We note that C is a continuum disconnecting C into 
two components, see Proposition 2.6 . Also note that the bounded component of 
C\C equals ft =^" 1 (Q). 

Next consider an increasing sequence Qi <s Q2 <e • • ■ <e ... Q of open squares, 
Qn = KQ, < Ai < A 2 • • • ->■ 1, so U^Li Qn = Q- Hereafter, the notation AQ 
stands for a square with the same center as Q but A-times smaller than Q. We will 
be dealing with the induced cells Oi (s f^2 <s • • • (s . . . O , where 0„ = h^ 1 (Q n ) , so 

U^Li ^« = ^- Choose one of the boundaries dQ n , n — 2, 3, This is a Jordan 

curve which separates Q„_! from dQ n+ \. By virtue of monotonicity of h: C C 
we see that C n := h~ 1 (dQ n ) is a continuum in 0„ +1 . This continuum separates 
dfin+i from fl n -i- Precisely, C n C Q n +i and the bounded component of C \ C n 
contains n„_i, because no point in 0„ = h^ 1 (Q n ) can be path-connected with 00 
within the open set C \ C n . 

Now supposed we are given a sequence {hj}°^ 1 of continuous monotone mappings 
hj : C — > C , such that h 3 ■ : C — > C , and /i(oo) = 00 . We assume that this sequence 
converges uniformly to /i: C — >• C. Denote the preimages of 9Q„ under /ij by = 
hJ 1 (dQ n ) , j = 1, 2, ... . These are continua disconnecting C into two components. 
The sets 

(3.1) ni:=hj\Q n ) 

are the bounded components of C\ C J n . Just the uniform convergence hj h yields 

lim sup dist(x, C n ) = . 
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Therefore, for sufficiently large j , say j j n , the continua C J n separate fl n -i from 
dfl n+ i. We record this fact as 

(3.2) H n _i C W n C n 3 n C n„+i , n ^ 2 and j ^ j„ . 

3.1.1. Cutting a cell. From now on the continuous monotone mapping h: C — > C 
will be the one that is obtained from the extension of h : X Y between circular 
domains, as in §2.6. Fix one of the boundary circles X\, . . . ,Xt C <9X, say Xi for 
some i = 1, . . . ,£, and recall that Ti = h(Xi) is the corresponding boundary circle 
in <9Y. We also recall the notation Xi for a reflection of X about Xi, and Yi for a 
reflection of Y about T^ Now consider an open square Q that intersects Yj along 
an open arc 7 = T, n Q. We assume that Q is small enough so it lies entirely in 
YU Ti U Yj. The arc 7 is a cross-cut of Q; it cuts Q into two connected subdomains 
Q n Y and Qn Y ; . We shall now discuss in some detail an analogous cross-cut of 
the cell f2 = /i _1 (Q). To this end we notice that the boundary map h: Xi ^> Tj 
defines an open subarc of Xi which lies in Q 

^{i£lj: h(x) G 7} . 

The endpoints of /3 belong to <9S1. Indeed, if x is an endpoint of /3 then there are 
points x v € j3 C which converge to x £ (3, thus h{x v ) £ 7. Passing to the limit 
we obtain h(x) € 7 \ 7 C dQ. Therefore, x £ h^ 1 (Q) = 0, meaning that x <G d£l. 

This is surely a geometric folklore that an arc j3 C (in a simply connected 
domain fi) whose endpoints lie in dtt splits f2 into two connected subdomains; 
namely, 

(3.3) f/=Xnfi and Xi n 
This latter subdomain will be of no interest to us. 

3.2. Hurwitz's theorem for p-harmonic mappings. Let us begin withg the 
classical Hurwitz theorem: 

Theorem 3.1. If a sequence of holomorphic functions f n : fl — > C converges c- 
uniformly to a holomorhic function f : O — > C, and f n (z) 7^ for all z G f2 and 
n = 1,2, ... , t/ien / is either equal identically to zero or f(z) ^ for all z G O. 

As a particular consequence of Hurwitz's theorem one obtains 

Theorem 3.2. If a sequence of conformal mappings ip n : Q — > C converges c- 
uniformly to ip: ft — >• C, i/ien (/? is either constant or conformal. 

In our arguments we will explore an adaptation of these two results to quasireg- 
ular mappings. It is immediate from Stoilow factorization of quasiregular map- 
pings [1, 3] that 

Corollary 3.3. In Theorem 3.1 holomorphic functions can be replaced by K- 
quasiregular mappings while in Theorem 3.2 conformal mappings can be replaced 
by K-quasiconformal homeomorphisms. 

In this section we aim to prove the following analogue of Hurwitz's theorem in 
a p-harmonic setting. 

Theorem 3.4. If a sequence of p-harmonic orientation-preserving homeomor- 
phisms ip n : Q — > C converges c-uniformly to (p: Q — > C, then either ip is a p- 
harmonic homeomorphism (actually ff 00 -diffeomorphism) or J v (z) = in O. 
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Proof. For basic properties of quasircgular mappings in relation to p -harmonic 
functions we refer the reader to [17] and for further reading to [22]. Recall from [33] 
that every p-harmonic function u: Cl — > R actually belongs to ^'"(Ct), where the 
largest integer k ^ 1 and the Holder exponent a € (0, 1] are determined through 
the equation 

7p - 6 + vV + 12p - 12 

fc + a = ^ . 

op — 6 

Thus, regardless of the exponent p, we have u £ ^ '°(fi) with a > 1/3. The 
complex gradient 

(3.4) / :=u z = ^ (u x - iu y ) , z = x + iy 



is a if-quasiregular mapping 



p-1 if p > 2 

(p-1)- 1 ifl<p<2. 

Actually, / satisfies the quasilinear elliptic equation 

df = 2-pH df f W 
{ ' ' dz 2p [f dz f dz 

In particular, every sequence {u n }^ =1 of p-harmonic functions converging c-uniformly 
on Cl, actually converges in ^^'"(Cl) . Then, by Corollary 3.3 , if the complex gradi- 
ents /" = u z are nowhere vanishing in Cl then / = u z is either identically zero or is 
also a nowhere vanishing function. Moreover, /" — > f uniformly on compact sub- 
sets together with all derivatives. Now consider the p-harmonic homeomorphisms 

ip n := u n +iv n : Cl^C 

and their limit 

(p:=u + iv ip:Cl—>C. 
Theorem 2.16 tells us that ip n are diffcomorphisms of positive Jacobian determinant, 
J (fin (z) — u™Vy — UyV™ > . Therefore, the complex gradients, 

/"(*):=< and g n {z):=v n z 

do not vanish in Cl . In view of Corollary 3.3 the limit functions f(z) — u z and 
g(z) = v z do not vanish as well, unless they are identically equal to zero. But 
in this latter case we would have J v = u x v y — u y v x = 0, proving Theorem 3.4. 
Therefore, let us assume that neither / =/= nor g =/= , everywhere. We aim now 
to show that also J v (z) ^ everywhere in CI. Let us notice in advance that once 
the inequality J ip (z) ^ is established, the map <p will be a local diffeomorphism. 
Even more, since this map is a c-uniform limit of homeomorphisms, it must be a 
global homeomorphism, by elementary topological considerations. 
Striving for a contradiction, assume that J v (z ) = at some point z Q e Cl, or 
equivalently, 

af{z Q ) + (3g(z ) = for some real numbers a, f3 ^ 0. 
We consider the complex functions 

F n (z) = af n (z) + Pg n {z) 

and their limit 

F(z)=af(z) + (3g(z), F(z o ) = 0. 
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Given any real coefficients a, f3 7^ we shall construct elliptic first order system of 
partial differential equations for F" . These equations will be of aid in determining 
that F(z) 7^ in 0, contradicting the equality F(zq) = 0. The derivation of the 
equation goes as follows. 

Using (3.5) and the analogous equation for g we compute 



F? = 



+ 



2-p 
2p 

2-p 
2p 

2-p 
2p 



f n df n f n df 71 



fn Q z fn Q z 



a + 



P 



2p 



dg n g n dg 
dz g 



dz 



g n Qpn g n Qpr, 

g n dz 
7™ _ 

fn g n 



r dz 

Qfn 

dz 



+ 



dz 



a 



Here we have 



and hence 



■11 



(af n + (3g n )f" - (a/« + pg»)f» 
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Thus we have the following first order (elliptic) inequality 



dF 7 



dz 



1 



dF n 



dz 



1 



2H 

P\\g r - 



This inequality can be viewed, equivalently, as a linear equation 
dF n 



(3.6) 



H n {z)^- +A n (z)F n , a.e.infl. 



dz dz 
with complex measurable coefficients satisfying: 

\fj, n (z)\ «; fc : 
|A n (z)| 2fc 





2 






1 




< 1 




P 




a 




1 






~P 




5™ 




dz 



e^c(tt). 



To see this latter inclusion we fix a compactly contained subdomain Ggfi. Recall 
that the functions g n (z) and their c-uniform limit g do not vanish. Thus we 
have a uniform bound from below |<7™(z)| m , for z e G, m being independent 
of n = 1,2,.... On the other hand the continuous functions ^-(2) converge 



uniformly on G , so we also have the uniform bound 



2kM\a\ 



or 

dz 



< M . This yields 



in G , as desired. 



Next we solve (uniquely) the nonhomogeneous Belrami equation 

(3.7) A" = [/i"(z)A™ — A n (z) ] Xg( z ) for a complex function A e^ 1 ' 5 ^ 2 ) 

where s > 2. We search the unknown function A"(z) in the form of the Cauchy 
transform of a complex density function w n € Jf s (M 2 ), with suppw" C G , namely 



""(f) df 
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The density function u>" will be found (uniquely) by solving a singular integral 
equation 

Here 5: Jz? s (C) if s (C) is the familiar Beurling-Ahlfors operator. Thus 

e> n = -(I-n Xa Sy'A^a 

We note that the operator 

is invertible for all 2 < s < 1 + j: , see [5]. We see that sup„||u; n ||^fs(c) < oo and 
infer that the family {X n }^ =1 is equicontinuous. Indeed, we have a uniform Holder 
estimate, 

\X n ( Zl ) - A"(z 2 )| < C s \ Zl - z 2 \ 1 -i\\cj n \\^ (C) . 

We may, and do assume, by passing to a subsequence if necessary, that X n (z) =4 A(z) 
uniformly in C. Therefore, the functions 

e r(2) F n (z) z4 e A(z) F(z) uniformly on G. 

On the other hand the nowhere vanishing functions H n (z) = e x F n satisfy a ho- 
mogeneous Beltrami equation 

fl£ - in G 

which is straightforward from (3.6) and (3.7). In other words H n are i-T-quasiregular 
mappings with K — y^f • Using Corollary 3.3 we see that the limit function H(z) — 
e x ( z ^F(z) docs not vanish anywhere in G, unless F = . This later case would 
mean that F(z) = af(z) + (3 g(z) = so, J v (z) = 0, which we have already ruled 
out. □ 

3.3. p-harmonic replacements. Here X and Y are circular domains. Let h <E 
Jf?p(K, Y), p > 2, and Q be an open square in Y + such that Q n Y is convex. Then 
the set 

U = r^QnY) = {x g X: fc(a;) eQflY} 

is a simply connected domain in X. We call it a cell. More specifically, U will be 
called an inner cell if Q C Y and a boundary cell if Q D dY ^ . Note that Q can 
intersect only the outer boundary circle of Y, since otherwise QnY would not be 
convex. 

Proposition 3.5. Let h <G Jf%,(K,Y) and U C X be a cell. Then there exists 
h* : X ^ Y , h* = h: <9X ^ dY, such that 

(i) h* e J^(X,Y) 
fit) /i* = ft: X \ £/ °-^> Y \ (Q n Y) 
(Hi) h* : U QC\Y is a p-harmonic diffeomorphism 

Moreover, 

(iv) ff |vfcT< // IW 
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Proof. We begin with a general consideration, without making explicit distinction 
between the two types of cells. Subsequently, additional details will be required for 
the case of a boundary cell. But at the end the proof will go in the same way for 
both cases. Recall that h is a weak limit of homeomorphisms hj € Jf? p (K,Y) . 

We conveniently extend each hj and h, as in §2.6, to obtain the monotone map- 
pings h j : C C and their limit h: C C. Note that {hj}JL 1 are homeo- 
morphisms in X as well as in the reflected domains X, , i = 1,2, ... ,n, but not 
necessarily along the boundary circles, so the mappings 

hj : X+ = X U Xi U • • • U X„ ^ Y U Yi U • • • U Y„ = Y+ 

are only continuous and monotone. Let 3Ej and Tj , for some i = 1, . . . ,n, be the 
outer components of <9X and dY, respectively. Thus Q C Y U Y, U Tj, which may 
or may not intersect Tj . Recall from §3.1 a sequence of open squares Q„ = A„Q, 
where < Ai < A2 <■■■—> 1, and the cells 

O n =^ 1 (Q„) CXUX.UI,, Oi (1 Q 2 (1 ••• (1 ...n = h- 1 (Q). 

Also recall the continua C n — hj (dQ n ), the bounded components 0^ = hj (Q n ) 
of C \ C n and the inclusions 

0„_i CflJC^C f2 n+ i , n ^ 2 for j > j n . 

In case of the inner cell U = h^ 1 (Q), the square Q lies entirely in Y, so the continua 
C n are Jordan curves, because hj — hj : fl^ ^> Q„ are homeomorphisms. Let us 
fix j, say j = j n , and introduce a sequence of homeomorphisms hj n : fi^™ Q„, 
n = 2, 3, . . . We shall return to these homeomorphisms after thorough discussion 
of a similar construction for the boundary cells. Let 

U = r'fQnf) = {x e X: eQnf} = onx. 

As observed in §3.1.1 this simply connected domain has resulted by cutting fl — 
h^ 1 (Q) along an arc in the boundary circle X = 3C,. This arc splits each of the cells 
f?2, . . . into the pair of two simply connected subdomains; the one which lies in 
X will remain of interest to us. We denote it by U n — f2„nX. Clearly, U^Lj U n = U. 
On the other hand, for j ^ j n we have simply connected domains fl 3 n — hJ 1 (Q n ) 
which are the inner complements of the continua C n — hj (dQ n ) C ^ n +i- But we 
will be interested only in a portion of fl n that is compactly contained in X. To this 
effect, we shear the set Q n Y along the circular part of its boundary. Precisely, 
for S > sufficiently small we consider the region Q n n Y s C Q„ n Y, where 
Y s = {y e Y: dist(y, dY) > 5}. The portion of fl J n that we are going to consider 
is a Jordan domain fl^ 5 = hj 1 (Q n n Y s ) C X with j = j n . Here 5 = 6 n > is 
chosen as follows. For 5 = we have 

nil = hj 1 (Q n nY) = njnx, j - j„. 

Since ftj : X Y is a homeomorphism, there is 5 = 5 n > such that the set 
Q£<5 = /iT 1 (Q n n Y" 5 ) contains Q,{ n X e ™, where X c ™ = {1 6 X: dist(:r, <9X) > e„}, 
say £„ = i . Let us write it as 

Jlj 5 - Z^ 1 (Q„ n Y 5 ) D 0„_! n X e " , where j = j n , 5 = S n and e„ = - . 
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Thus we have a homeomorphism 

(3.8) h jn : f!^ 5 " ^Q„nY 5 " 

of a Jordan domain 17^™ ' Sn onto the convex domain <Q>„ n Y Sn . 
From now on the remaining part of the proof treats both types of cells in the same 
way. Simply, in case of the inner cells we have il 3 r ^- Sn — and <Q>„ n Y Sn = <Q>„. 
We now appeal to Theorem 2.16. Accordingly, each homeomorphism in (3.8) can 
be replaced by a p-harmonic diffeomorphism, denoted by 

~ h = (h jn :X\ fl^Sn _> Y \ (Q„ n Y Sn ) 

1 p- harmonic diffeomorphism of fl^' Sn ^^^^ Q„ n Y Sn . 

By Theorem 2.15 we still have hj n £ Jf? p (X,Y) and 

£ x\hj n ] < £x[hj n ] , for j 1 < j 2 < > oo . 

Now the desired mapping h* £ Jtfp(X,Y) is defined to be the weak limit of hj n . 
Recall that homeomorphisms in ^,(X,Y) admit a continuous extension to the 
closure of X , with values in Y . Furthermore, hj n converges uniformly to h* on 
X , because of a uniform bound of energy and Lemma 2.10. For x £ X \ U we have 
hj n = hj n — > h(x), so h* = h on X \ U. Moreover, h* is a p-harmonic mapping on 
U = ft n X. By Theorem 2.15 we sec that h* £ h + W hp {U) and 

£ x [h*} < £ x [h}. 

It remains to show that h* : U n X Q n Y is a diffeomorphism. To this effect, 
let U' be any compactly contained subdomain of U. We note that U' C Q J ^' Sn 
for sufficiently large n. By Hurwitz type Theorem 3.1 h* : U' — > C is either a 
diffeomorphism or its Jacobian determinant vanishes identically on U'. However, 
this latter case is easily ruled out by computing the integral of the Jacobian over 
U. For this purpose we appeal to the equation (2.13), 

/ J(x,h*)dx = lim / J(x,hj ) dx > lim / J(x,hj )dx 

= lim |Q„n Y 5 "| = |Qn Y| > 0. 

In particular, j u , J(x, h*) dx > for sufficiently large subdomains U' <s U. In con- 
clusion, h* : U' — > C is a diffeomorphism. Letting U' approach U we conclude that 
h* : U —¥ C is a diffeomorphism. Recall that homeomorphisms take O^"'' 5 ' 1 C {/ 
onto Q„nY d " C Q„nY and converge uniformly on U to h* . Hence, by elementary 
topological arguments we see that h* : U Q n Y is a diffeomorphism. 

3.4. Further replacements. A key step in constructing homeomorphisms h* : 
X Y for Theorem 1.4 is the following 

Proposition 3.6. Let X and Y be circular domains bounded by circles X\, . . . ,Xg 
and T 1; . . . , respectively, and h £ <^,(X, Y). Select a pair of boundary circles 
X = Xi, T = Tj and recaH i/iai /i: Xj Tj . Consider the open region in X 

O = r^YUT) = {x £ X: /i(ar) e Y U T} . 

Then for every e > t/iere exists H : X — >• Y \ T in Jf p (H,Y) which satisfies the 
following conditions: 
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(a) H = h : dX °-^> dY , also H = h: X \ Q -> dY \ T 

(b) H : VL Y is a homeomorphism 

(c) \\H - h\\sg P (X) < e 

Of course 7? depends on e. However, to ease the writing, we suppress the explicit 
dependence on e in the notation of H but return to it later on. 

Proof. We shall make use of Proposition 3.5, which requires that T be the outer 
boundary of Y . Let us first demonstrate how Proposition 3.6 can be reduced to 
this case. 

Step 1 (Reduction to the case of outer boundary) Given a pair (£j,Yj) of 
boundary components, we recall that h: X, Yj . We explore the reflections 
ifi : X X, and ipi : Y Yj from §2.6. In this way X, and Yj become 
outer boundaries of the circular domains Xj and Yj , respectively. Consider the 
mapping 

V = i> i oho<p i :X i ^Y i , tfe J^(Xi,Yi). 
Assuming that Proposition 3.6 holds in case of outer boundaries we infer that given 
any e l > (to be revealed later on) there exists H l : X, -> Yj \ T l , W e J^(Xj, Yj) 
such that 

• W = h l : % \ Q! -> dY, \ T 4 , where W = {x e X 4 : K l {x) e¥,U TJ 

• H l : 51* OJ ^r Yj is a homeomorphism 

• \\W - WW&pQLi) < ^ 

We then find the desired mapping H e Jf p (X,Y) by setting 

= tpi o H i o ^ : X Y 

Indeed, the conditions (a) and (b) are obviously satisfied. Regarding the condi- 
tion (c) we note that ipi and ipi are ^°°-smooth diffcomorphism so 

\\H - /i|Up(x) = ||V>i ° H% o ^ - Vi o h l o ^ill^p(x) ^ C\\H' 1 - h'WsgpCKt) < Ce l = e 

where C depends only on the Lipschitz constants of the reflections; thus we take 
e* equal to e/C. The interested reader may wish to observe that this transition 
(by reflecting about circles) retains control over \\H — h\\gg P /x), but at the sacrifice 
of loosing the sharp inequality ||i?||^P(x) ^ H^H^p(x)- However, we shall have no 
need of this sharp inequality. 

Step 2 (Three families of squares) 

Lemma 3.7. Let ¥ be a compact subset in C and p > 0. There exist three families 
A = {^4q}2Li, B = {-BsjjgLi and C = {C 7 }^L 1; each of them consists of disjoint 
open squares in C \ F of diameter p or less, such that 

oo oo oo 

(J A a U [jBpU (J C 7 =C\F. 

a = l 0=1 7=1 

The construction of such families presents no difficulty We shall not bother the 
reader with an explicit construction of A,B and C. In our application the set F 
will be the union of all bounded components of C \ Y. We will be dealing only with 
those open squares Q which intersect Y, so Q n Y will be nonempty and 

oo oo oo 

(3.9) |J^ q u|J%u|Jc 7 dyut. 

a=l 0=1 7=1 
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To each family A, B and C there will correspond a mesh of cells in X . We begin 
with the family A . 

Step 3 (p-harmonic replacement in A-cells) Fix p > for the moment, small 
enough so that if an open square Q of diameter p or less intersects the outer bound- 
ary T = Tj C <9Y then it lies entirely in Y U Y, U Yj. Hence, in particular, Qfl¥ 
is convex. Consider the mesh of open cells in X 

U a = h~ 1 (A a n Y) = {x G X: h(x) G A a n Y} . 

It should be noted that U a is empty whenever A a n Y = 0. 




FIGURE 1. A mesh of cells in X and squares in Y 



We now appeal to Proposition 3.5 which provides us with a mapping that we 
denote by h A ■ X Y , such that 

(i) h A G J%,(X, Y) , h A = h:dX ^ «9Y 

OO OO 

(ii) h A = h:X\ |J C/ Q °-^Y\ |J(A Q nY) 

a— 1 a— 1 

oo oo 

(iii) hj± : M £7 Q (A a n Y) is a p-harmonic diffcomorphism. 

a— 1 a— 1 

(iv) JJWh A \ p < jjW- 

The actual construction of h A is accomplished step by step. First, using Proposi- 
tion 3.5 for h and U = U\, wc obtain what has been denoted by h*. In the second 
step we use Proposition 3.5 for h* in place of h and for U = U2 to obtain h**. 
This process continues indefinitely resulting in a sequence of mappings converging 
weakly to h A . Conditions (ii) (iv) are clearly satisfied. However, the condition (i) 
calls for some explanation. Upon every finite number of steps we obtain a mapping 
h**-* G J^(X,Y). These are weak ^'Mimits of homeomorphism in Jf? p (X,Y). 
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On the other hand, as the number of steps approaches oo, the mappings h**'"* 
converge weakly in W 1,P (K, Y) to h A . Hence h A e Jf p (X,Y) , as desired. 

It should be recalled at this point that the family A of squares actually depends 
on p > . We write it as A = A p . Therefore, for each p we have mappings h Ap e 
Jf p (X,Y) which satisfy conditions (ii)- (iv). Let us now take a close look at the 
mappings h A when p approaches 0. First, given x € X, we see that h A (x) = h{x) 
if x does not belong to any of the cells U a . If, however, x € U a for some a then 

h A (x) E h(U a ) = A a n Y . 

In either case \h a (x) — h(x)\ ^ p. We then infer that \\h Ap — h\\<#(x) ^ P, which 
means that h Ae =4 h uniformly on X as p — > 0. On the other hand the energy of 
h As is not greater than that of h. Thus h Ag — > h, weakly in W 1,v {%). Hence 

£x[h] < lim £x[h Ap ] sC lim £ x [h] = £x[h] . 
p— > p— >o 

This shows that in fact lim £x\h A ) = £x\h]. Consequently, h A , — > /i strongly in 

p— > o 

M P {X). In summary, given any e > 0, we can choose and fix p > for the family 
A = {Aa}^ to satisfy 

(3.10) \\h A - h\y p{X) s? * e 

Step 4 (p -harmonic replacements in B -cells) We repeat the procedure in Step 
3, with /i^i in place of h and with the family B in place of A; with the same e > 0. 
Thus, we consider the mesh of open cells in X 

V fj = h A 1 (B f3 nY) = {xeX: h A (x)eB nY}, 

and we obtain a mapping, denoted by h AB : X Y, such that 

(i) h AB G ^(X,Y), ^ B = h A = h: <9Y 

OO OO 

(ii) /MS = /U:X\ [J^ A¥\ (J (Bp n Y) 

0=1 0=1 

oo oo 

(hi) /i^g : ^ Vg (Bp n Y) is a p-harmonic diffeomorphism. 

0=1 

(iv) / f \vh AB \p < // iv^r < // iv^r 

././x ./ix ././x 

The inequality parallel to (3.10) reads as 

\\h A B - h A \\& P (x) ^ 3 e • 

Step 5 (p-harmonic replacements in C -cells) The above procedure of p- harmonic 
replacements is finally carried over to the mappings h A B ■ We consider the mesh of 
open cells in X, 

W 1 = h^ B (C^ n Y) = {x e X: h AB (x) e C* 7 n Y} . 
The replacements result in the mapping, denoted by h AB c '■ X a -^>- Y, such that 

(i) h ABC G ^(X, Y) , h ABC = h AB = h A = h: 9X ^ OY 

OO OO 

(ii) h ABC = h AB :X\ |J W 1 °-^Y\ \J(C y nY) 

7=1 0=1 
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oo oo 

(iii) h A sc'- [J W 1 [J (C 7 n Y) is a p-harmonic diffeomorphism. 

7=1 fi=l 

(iv) f [ \Vh ABC \ p < ff \Vh AB \ p < ff \Vhf 
Moreover, 

(3-11) \\fiABC - hAB\\sep(x) < 2 e - 

Step 6 f 77ie desired mapping H = h A bc) Upon three consecutive replacement 
procedures in the *4-cells, S-cells and C-cells we finally arrive at the mapping h A Bc- 
This map turns out to satisfy all the conditions asserted in Proposition 3.6, so we 
denote it by 

Let us verify those conditions. Condition (c) is obvious; 

\\H - ft 1 1 a* (X) ^ \\hABC - ft^e||«p(x) + \\hAB - ft^||«p(x) + ||ft^t - ft||«p(x) 
e e e 

^3 + 3 + 3= e - 

Verification of other conditions involves elementary set theoretical considerations. 
Let use include them for completeness; some do not seem immediate at all. Let us 
begin with the following identity 

oo oo oo 

(3.12) U u a u U vpu |J w 7 = n 

a=l 0=1 7=1 

where ft = {x G X: ft(z) G YU T}. 

Proof of (3.12). Take any point in X such that 

oo oo oo 

X |J U a U |J Vf, U |J TU 7 . 

a=l 0=1 7=1 

Since x G" U^Li ^a, we see that h(x) = h A (x). Since x £ UsLi we see that 
ft_4(a:) = h AB {x). Since a; ^ U^Li ^7' we see tnat = h AB c{x) = H(x). 

Hence y := ft(a;) = ft .4(2;) = H A b{x) = h A Bc{x)- Moreover, we have 

oo oo oo 

y = h(x) £ |J (A Q nY) , y = h A (x) £ {J (B^nY) and y = /^(x) £ (J (C 7 nY) . 

a=l 0=1 7=1 

Thus 

(OO OO oo \ 

|J A a U [J BpU |Jc 7 nf. 
a=l /3=1 7=1 y 

On the other hand, by the formula (3.9) , h(x) ^YUTsoi^J]. 

Conversely, suppose x G U^Li U U/jLi ^8 u U^i W-y ■ Then we have three 
cases 

Case lie |XLi C 7 "' which yields h(x) € A Q n Y C Y U T, so x G ft. 
Case 2 x ^ \J^ =1 U a and x G U/3L1 Xs- This yields 

h(x) = h A (x) G Bp n Y G YU T, ie!!. 
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Case 3 x £ U^Li U a and x £ UJ3L1 Vp , so x e U^Li This yields 

/i(i)=Mi) = /ii B (i)eC 7 nfcYuT, a;eO. 

completing the proof of (3.12). □ 

A care needs to be excercised in showing that H : Q — > Y is injective. For this, 
consider two distinct points 

oo oo oo 

xi,x 2 e n = |J t/a U |J ^ U |J w 1 

a = l 0=1 7=1 

and suppose, to the contrary, that y := H(x 1 ) = H(x 2 ). Since the mapping 
H = h A BC : U"=i W 7 U~=i (C 7 n Y) is injective and 

OO oo 

h A BC = h AB : X\ |J IT 7 ^Y\ (J(C 7 nY) 

7=1 7=1 

it follows that xi,x 2 ^ Ut=i W-y- Thus we are reduced to the case xi,x 2 € 
UsLi ^ UsLi ^8 anc ^ 2/ = ^Ab( x i) = h AB (x 2 ). In exactly the same way we 
infer that Xi,x 2 $ U/3=i Vp- Thus we are further reduced to the case 

oo 

xi,x 2 e\JU a and y = h A (xi) = h A (x 2 ) . 

a=l 

Just as before Xi,x 2 U^Li U a - We arrive at the contradiction 



xi,x 2 # \J U a U \J VpU \J W 7 

a=l 0=1 7=1 

finishing the argument for injectivity. 

To see that H : f2 '-^ Y, we argue in similar way; namely, let 

oo oo oo 

a;efl=U[/ a u|J^u[Jlf 7 . 

a = l 0=1 7=1 

If xe U^Li W, n then 

OO 

y:=H(x)=h ABC E (J C 7 H¥c Y. 

7=1 

Thus suppose a; ^ U7I1 ^7> so = h A s{x) and x e U^Li ^« u U/3L1 Vp- ^ 

x e U/3L1 tnen ?/ — h AB e U/3L1 ^n¥cY. Then suppose that a; UJ3L1 ^81 
so x G U/3L1 U a - This yields y = h A (x) e U^Li 4nYcY which completes the 
proof of the inclusion H(£l) C Y. 

Lastly, to see surjectivity H: Y we recall that H = h A sc'- X Y and 

H = h: dX ^> dY. Hence ff(X) D Y. Take any y <= Y. There exists iel such 
that H(x) = y. Suppose that x £ = U~=i ^ U (j^=i ^ U U^Li W7. This leads 
to a clear contradiction, 

y = H(x) = h ABC = h AB = h A = h(x)€Y, and x e /i _1 (Y U T) = SI . 

The proof of Proposition 3.6 is complete. □ 
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3.5. Proof of Theorem 1.4 for multiple connected domains. Given h G 
Jf p (K, Y) , we define by induction a chain of mappings H , Hi , . . . , Hi G J4? P (X, Y) , 
setting H = h. Suppose that H k _i G J^ p (X,Y) is already defined for some 
1 ^ k ^ i. To define Hk we consider the boundary circles X k C <9X and T k C <9Y, 
and recall that h: Xk T k . We consider the open regions 

n k = H-\(Y UT k ) = {xeX: H k _i{x) G Y U T k } . 

By Propposition 3.6 applied to H k -i in place of h we obtain a mapping H , which 
for our purpose will be denoted by H k G Jf° p (X,Y). We have 

(i) H k =H k _ 1 :X\il k ->dY\T k , H k = H k _ x : dX ^ 8Y 

(ii) H k : Q k Y is a homeomorphism 
(hi) \\H k - H k _i\\^ P (x) < e 

The desired sequence ft-* G Jt? p (X, Y) of homeomorphisms /i* : X Y converging 
to h strongly in #' 1,P (X, Y) will be defined as follows. Let a sequence of positive 
numbers Cj — > be chosen and fixed. Then we set h* = Hi, where Hi is the last 
link in the chain H , Hi, . . . , Hi, with e = ej. We need to verify four claims. 
Claim 1 We have 

fiiUfi 2 U'"Ufi£ = X. 

Here each Cl k is contained in X, so we need only show that every i£X lies in 
one of those domains. Suppose, to the contrary, that x Oi U • • • Ufi^. Since x G" Q,i 
we see that Hi(x) = Hi_i(x), and by induction that Hi_i(x) — Hi- 2 {x) — ■ ■ ■ = 
H (x) = h(x). These equalities, together with the definition of Cli, . . . , Slf yield: 

h(x)=H (x) ^ YUTi 
h{x)=Hi(x) £ YUT 2 

h(x) = H t -\{x) ^¥UT< 

Thus /i(i)^¥UTiU'"UT,=Y,a clear contradiction. 

Claim 2 H e :X°-^Y. We have H e (X) D H e (Sl e ) D Y. Now take x G X 
and assume, to the contrary, that Hi(x) G dY. Since Hi : Qi Y it follows that 
x £ fli which in turn means that Hi(x) = Ht_i{x) G dY. But : f2^_i Y, 
so x £ fli-i. Continuing in this way we conclude that x Cli. This summarizes as 
x £ Qi U • • • U fli = X, a clear contradiction. 

Claim 3 H f : X Y is injective. Suppose, to the contrary, that Hi(xi) = 
Hi{x2) = y G Y for some xi 7^ £2. Since -ff^: Y is a homeomorphism it is 

impossible that both x\ and x 2 belong to Qi. It is also impossible that xi G £li and 
x 2 ^ flf, since otherwise we would have Hg(xi) G Y and Hi(x 2 ) G <9Y \ T^. This 
leaves us with the only possibility that xi, x 2 ^ Ui, so y G <9Y \ T^, contradicting 
Claim 2. 

Claim 4 Obviously, we have 

e 

\\tij - h\\®p{X) < y^ll-^fc - #fc-i||««>(x) < -> 0. 

fe=0 



□ 
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4. The case of simply connected domains 

Theorem 1.4 fails when t = 1 and p = 2. To see this consider a sequence of 
Mobius transformations h k ■ D D of the unit disk onto itself, 

(4.1) h k (z) = ±±*L, afc =^-H, h k (l) = l, h k (-l) = -l. 

1 + CLkZ K 

All the mappings have the same finite Dirichlet energy. 

£n[h k ] = f[\Dh k \ 2 = 2 [f J(z,h k )dz = 2w. 
■I .1 : J J : 

They converge to the constant mapping h{z) = 1, c-uniformly and weakly in 
W 1,2 {&). Obviously h(z) cannot be obtained as a strong limit in W 1,2 {p) of homc- 
omorphisms of D onto itself. Nevertheless, continuous monotone extensions of h k 
outside D are still available, but at the sacrifice of loosing uniform bounds (inde- 
pendent of k) of the Dirichlet energy in any domain QdD. Indeed, otherwise we 
would have a uniform bound for the modulus of continuity, which in turn would 
imply that the boundary mappings h k : dH 9D converge uniformly; this is not 
the case. Thus we see that reflecting D about the boundary circle cannot help. 

4.1. Proof of Theorem 1.4 for simply connected domains, p > 2. When 
the Sobolev exponent is greater than 2 the homeomorphisms h k <G J^,(X, Y) and 
their weak limit h € J4? P (K,Y) enjoy uniform modulus of continuity, see Lemma 
2.10. The reduction to the case of unit disks, X = D and Y = D, can be made 
via a bi-Lipschitz transformation as in Section §2.5. This reduction is needed only 
to make the target domain convex, so as to apply the p-harmonic replacements. 
The point is that there is no need to reflect h k or h about the boundaries of X and 
Y; radial extensions outside the disks (for topological purposes) work just as well. 
The rest of the proof of Theorem 1.4 runs as in the case I ^ 2 , with hardly any 
changes. 

4.2. Simply connected domains with a puncture. Let us now take on stage 
homeomorphisms h k : X Y, h k e ^(X, Y) , between simply connected Lips- 
chitz domains. Assume that the mappings h k are fixed at a point x Q e X, namely. 

h(x ) — y , k — 1,2, ... , where x G X and j e¥. 

Because of this normalization we have an extension h k : C "-^^ C without loos- 
ing uniform bounds in a neighborhood of X. Indeed, with the aid of bi-Lipschitz 
transformation, see Section §2.5 , we are reduced to the mappings h k : D D 
between unit disks such that h k (0) = 0. With this normalization the mappings 
can be interpreted as homeomorphisms h k : D D between punctured disks 
D = D \ {0}. By Theorem 8.1 in [35] we infer that 

(4.2) dist 2 [/*(*), 3D] < ° ' £ ° [hk] 



log {, 



e + 



dist(z,dO) 



where C is a universal constant. Now consider the extension h: C — > C, by the 
formula 

-l 

for \z\ > 1 . 



h k {z) = h k (l/z) 
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We have £o[/ifc] ^ M for some constant M < oo and all k = 1, 2, . . . . Choose and 
fix r > 1 sufficiently close to 1 so that 



C-M 1 

^4 



log(e+^) 
Inequality (4.2) reads as 

" + Tift) lo g( e +^l) 4 

Hence | h k (l/z) | > 1/2 , which gives us the desired uniform bound of the Dirichlet 
energy 

\Vh k \r^ C- II \Wh k \P, k = l,2, 



I In, 



The remaining arguments for the proof of Theorem 1.4 are the same as in the mul- 
tiply connected case. Let us summarize these findings by the following statement, 

Theorem 4.1. Let X and Y be simply connected Lipschitz domains and leth k : X 

Y be homeomorphisms in W 1,2 (H,Y) normalized by h k (x ) = y Q , for some x Q e X 
and y € Y, and all k = 1, 2, . . . Suppose h k converge weakly in # /1 ' 2 (X, Y) to a 
mapping h G #' 1 ' 2 (X, Y). ITien i/iere exists a sequence of -diffeomorphisms 

h* k :X°-^Y, h* k eh + W ^ 2 (X,Y) 

converging to h strongly in #' 1 ' 2 (X, Y). Moreover, 

(4.3) h* k (x )=y for all fc = 1,2,... 

Remark 4.2. The reader may be concerned that our proof of Theorem 1.4 does 
not guarantee that /i^(a; ) = y Q - However, since h k (x ) — > h{x ) = y , one can 
alter slightly each h* k , by a smooth change of variables in X, to achieve such a 
normalization. 

4.3. A normalization at three boundary points. Consider homeomorphisms 
h k : X — > Y between bounded simply connected Lipschitz domains in the Sobolev 
space #^ 1,2 (X, Y), normalized by three conditions at the boundary: h{x\) = y\, 
h{x2) — 2/2 and h(xs) = 2/3 , where x\,X2,Xz are given distinct points in <9X and 
2/1)2/2,2/3 € dY . We assume that both triples (xi, X2, #3) and (2/1, 2/2, 2/3) are posi- 
tively oriented along <9X and 9Y, respectively. We transform the domains X and 

Y via bi-Lipschitz mappings into equilateral triangles in which (x\, X2, £3) and 
(2/1)2/2,2/3) are their vertices. As in the multiply connected case, we will need to 
extend the mappings to a neighborhood X + D X . To build this neighborhood, we 
reflect X about its sides. This results in three adjacent triangles. We then reflect 
the adjacent triangles. Eventually we come up with 12 equilateral triangles sur- 
rounding X which together with X make the neighborhood X + . The same reflec- 
tion procedure is made for Y. We then extend each homeomorphism h k : X <i; ^> Y 
in accordance with the above reflections and define continuous monotone extensions 
h k : X + Y + . The way the extensions are made secures a uniform energy bound 



ff |v£ fe | 2 = 13 [[ |V/i fe | 2 
J Jx + J Jx 



< M , for some M and all k = 1, 2, . . . 
By the same arguments as for Theorem 1.4 we draw the following conclusion 
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Theorem 4.3. Let X and Y be simply connected Lipschitz domains and let hk ■ X — 
Y be homeomorphisms of Sobolev class # /1 ^ 2 (X,Y). 1 Suppose h k are normalized 
at three boundary points; 

h{x{) = yi, h(x 2 ) = y2, h(x 3 )=y 3 

and that hk converge weakly in W 1,2 (1L,Y) to a mapping h e #" 1,2 (X,Y). Then 
there exist ( ^' oc -diffeomorphisms 

h* k :X°-^Y, h* k e h + W Q h2 (X,Y) 

converging to h strongly in # /1 ' 2 (X, Y) and uniformly on X. 
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